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We consider the problem of the validity of a statistical mechanical description of two-site en- 
tanglement in an infinite spin chain described by the XY model Hamiltonian. We show that the 
two-site entanglement of the state, evolved from the initial equilibrium state, after a change of the 
magnetic field, does not approach its equilibrium value. This suggests that two-site entanglement, 
like (single-site) magnetization, is a nonergodic quantity in this model. Moreover we show that these 
two nonergodic quantities behave in a complementary way. 
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Entanglement plays a key role in the rapidly grow- 
ing field of quantum information processing |]| , where so 
far, most of the discussions of entanglement concern few 
body systems. Recently, however, the properties of en- 
tanglement has been also used to study and understand 
behavior of "complex" quantum systems. For instance, 
entanglement in quantum many body systems such as 
spin chains or Bose-Einstein condensates were studied 
(see e.g. |2j,|sJ,|^). The role of "entanglement length" in 
quantum phase transitions was pointed out in Ref. 

The aim of this paper is to deal with the statistical 
properties of entanglement, or rather the validity of a sta- 
tistical mechanical approach to entanglement in realistic 
many body systems. The relation of entanglement with 
important notions in statistical mechanics is important 
to understand the behavior of entanglement in quantum 
many body systems. We will restrict ourselves to one- 
dimensional infinite spin systems. Entanglement in such 
systems was considered e.g. in Ref s. 0LiE|- These pa- 
pers deal exclusively with the properties of entanglement 
of an infinite spin chain, which is either in a ground state 
or in a thermal equilibrium state. Such states are "static" 
states of the system, as they do not explicitly depend on 
time. Exceptions include the recent studies of entangle- 
ment in quantum dynamics, which allowed to formulate 
novel types of unprecedentedly efficient numerical codes 
for simulation of quantum systems Q (see also H,0). 

The validity of a statistical mechanical description of 
a quantity, characterizing a physical system, depends on 
the behavior of that quantity as the system evolves in 
time. More precisely, a necessary condition for validity of 
statistical mechanical description of a physical quantity 
is that it must be "ergodic" . A physical quantity is said 
to be ergodic if the time average of the quantity matches 
its ensemble average. This is usually a hard question to 
check, even for classical systems. However, an indication 
of whether a given quantity is ergodic or not, can be 
obtained by comparing the time evolved state with the 
equilibrium state. 

To deal with such questions, we will therefore consider 
real time evolution of an infinite one-dimensional spin 
chain. We will suppose that the spin chain is described 
by the so-called XY model (see Eq. (QJ below). We 
will show that for large times, the nearest neighbor en- 
tanglement of the evolved states, does not approach to 



the corresponding entanglement of the equilibrium state. 
This suggests that nearest neighbor entanglement in this 
spin system is nonergodic [Toj , indicating that one cannot 
describe entanglement in such models by equilibrium sta- 
tistical mechanics. Temporal dynamics of entanglement 
in spin systems is important in quantum information and 
computation tasks. For instance, the temporal dynamics 
of entanglement in a spin system has been used in Ref. 
to obtain the one-way quantum computer. For the 
infinite spin chain in the XY model, it is known that the 
sing le-site magnetization is also a nonergodic observable 
We will show that (two-site) entanglement has a 
complementary temporal behavior to (single-site) mag- 
netization in this model. We will also show that the 
(two-site) entanglement and magnetization of the evolved 
state saturates for low temperatures, for a given time. 

The XY model (see Eq. QJ) that we consider here is 
integrable 01 ■ For classical systems, integrable models 
possess a large number of constants of motion, and are 
usually not ergodic (see e.g. However, such ob- 

servations are valid for the case of systems that have a 
finite (and usually small) number of degrees of freedom. 
In our case, we deal with an infinite spin chain, and our 
considerations are distinctly quantum. Hence, questions 
about ergodicity, or its absence are nontrivial. 

A one dimensional spin system (a one-dimensional ar- 
ray (lattice) of spin- 1/2 particles) with nearest neigh- 
bor interactions is described by a (dimcnsionlcss) Hamil- 
tonian of the form H mt = ^(ASfSf^ + BSfSf +1 + 
CSfSf +1 ), where the Sf, Sf, Sf are one-half of the Pauli 
spin matrices o~ x , o~ v , o~ z at the i-th site of the array, 
and A, B, C are coupling constants. Here, we take A ^ B 
and C = 0. We introduce also an external field into the 
Hamiltonian, so that the total Hamiltonian of the system 
is H(t) = Hi n t — h{t)H mag . To ensure that this external 
field has nontrivial effects on the evolution, we must have 
[Hint, H ma g] / 0. The simplest way in which this can be 



effected is by choosing H„ 



Ei Sf, and A= 1+7, 



B = 1 — 7, 7 ^ 0. In this way we arrive at the XY model 
in the transverse field 



H(t) 



E 
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[(l + 7 )SfSf +1 + (l 



i)s?s* +1 



h(t)S? 



(1) 

where 7 =/= 0. In the following, we set h = 1. A spin 
chain whose dynamics is described by H(t), is said to 
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be described by the "XY model". For a finite number, 
N, of spins, we assume a periodic boundary condition, 
Sn+i — Si. Ultimately, we will be interested in the 
thermodynamic limit N — > oo. Such systems can be 
realized in atomic gas in an optical lattice (see e.g. [l4j)- 
At a given time t, the (thermal) equilibrium state 
is given by p e ^{t) = exp(—(3H(t))/Z, where Z = 
tr(eacp(-/3.ff(f))). Here (3 = 1/kT, where k is the Boltz- 
mann constant, and T denotes the (absolute) tempera- 
ture. To consider questions about ergodicity, we will be 
interested in the behavior of the evolved state. The evo- 
lution is governed by the Hamiltonian H(t), from a given 
initial state. Since we will compare the properties of the 
evolved state, with those of the equilibrium state, it is 
natural to suppose that the initial state (at t = 0), is the 
equilibrium state at t = 0. We denote the evolved state 
by p a (t), where the suffix corresponds to the tempera- 
ture of the initial equilibrium state p^(0) = p a {0)- For 
later times, properties of the evolved state p a (t) will be 
compared to those of the equilibrium state p e ^ (t) , so that 
they have the same energies: 

tv(H(t)p a (t)) = tv(H(t)p eq (t)). (2) 

For simplicity, we will consider the case of sudden switch 
of the field h(t), as h(t) = a, for t < 0, and = b, for t > 0. 

Both p a (t) and p C p(t) are states of an infinite num- 
ber of spin-f/2 particles. For our purposes, it will be 
sufficient to consider single-site and two-site density ma- 
trices. Let us first consider the single-site density matrix 
for the state p e ^(t). By symmetry, the single-site density 
matrices of the chain are all the same. We will denote 
it by pl q (t) (hiding the suffix /3). Now p eq *(t) = p\ q {t), 
when the complex conjugation is taken in the compu- 
tational basis, which (for each site) is the eigenbasis of 
Pauli matrix a z . Therefore tr(S v p^ 9 (t)) = 0. Moreover 
the Hamiltonian H{t) has the global phase flip symmetry 
([H, U t S z ] = 0), from which it follows that ti(S x p\ q {t)) = 
0. Consequently, the single-site density matrix of the 
equilibrium state is of the form p eq {t) = \l + 2M^{t)S\ 
where I is the 2x2 identity matrix. The evolved state 
does not necessarily have the property of being equal to 
its complex conjugation, and consideration of the global 
phase flip symmetry is complicated by the fact that the 
Hamiltonian is explicitly dependent on time. However, 
using the Wick's theorem, as in 0, 0, 0] , the single- 
site density of the evolved state turns out to be of the 
form pi(t) = \l + 2M z {t)S z . 

For the case of the two-site density matrix p^(t) of 
the equilibrium state p e p(t), we can again use the global 
phase flip symmetry and the fact that it is equal to its 
complex conjugate, so that it is of the form 

pll(t) = -I ® 1+ M e z q {t)(S z ®I + I®S Z ) 



E 



T* q (t)S j ®S j , (3) 



where the correlation functions, T^ q (t), are defined as 
T-?(t) = itv(S j <g> SipH(t)), j = x,y,z. In the case 



of the two-site density matrix of the evolved state, the 
yz and zx correlations are absent (via use of the Wick's 
theorem). However the xy correlations does not vanish, 
just as the xx, yy, and zz correlations. Thus the two-site 
density matrix of the evolved state p a (t) is of the form 



Pu(t) = jl ® I + M z {t){S z ® I + I ® S') 
+T xy {t)(S x g> S v + S y ® S x ) + J2 T n(t)S j ®5 J ',(4) 



where the correlation functions, Tjk(t), are defined as 
T jk (t) = 4tr(^ ® S k p 12 (t)), j, k = x,y, z. 

Let us consider the (single-site) magnetization and the 
(two-site) entanglement in an infinite chain for the two 
states of interest. The magnetization of the equilibrium 
state is M^(t) = tr(5 2 p\ q '(*)) , while that for the evolved 
state is M z (t) — tr(S z pi(t)). For studying entanglement 
(Hit Il7j , we will use logarithmic negativity (LN) 01 as 
our measure of entanglement. LN of a bipartite state 
Pab is defined as E N (p AB ) = log 2 where ||.||i is 



the trace norm, and p T j^ B denotes t he p artial transpose 
of pab with respect to the A-part [Ijj. Note that the 
two-site density matrices in our case acts on C 2 ® C 2 . 
Consequently, a positive value of the LN implies that the 
state is entangled and distillable Il9l|20l . while En = 
implies that the state is separable |19| . 

The (single-site) magnetizations of the equilibrium 
state and the evolved state are respectively given by 
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tanh(i/3A(/i(t))) 

' MM) 



(h(t)-cos<f>), (5) 



and 



M z {t) 



d4> 



tanh(i/3A(a)) 



2tt J " t A(a)A 2 (fe) 
x [cos(2A(6)t) 7 2 (a-6) sin d> 
— (cos0 — 6)[(cos0 — a)(cos(j) - 



b) + 7 2 sin 2 i 



(6) 



where A(a) and A(6) can be obtained from A(h(t)) = 
[7 2 sin 2 <j) + (h(t) - cos 0) 2 ] 5 . 

For definiteness, let us consider the case where 7 = 
0.5, a = 200, a = 0.5, b = 0. For these values of the 
parameters, Eq. J2J gives f3 w 3.9. The equilibrium 
magnetization M z q (t) has an initial value of « 0T48328 
and then jumps down to zero for all later times. The 
magnetization M z (t) of the evolved state, on the other 
hand, is an oscillating function (see Fig. QJ, which of 
course starts at the same initial value [21| as M z q {t), 
but remains positive for long times. This suggests that 
the magnetization of the spin system in the XY model is 
nonergodic [HI (cf. HI). 

The nearest ne ighb or correlations for the equilibrium 
state are given by (HlH T x e |(i) = -G eq (-1, i), T$(t) = 
~G eq (l,t), TH(t) = 4[M z (t)} 2 - G e o(l,t)G e *(-l,t), 



where G eq (R,t), for R — ±1, are given by 



G eq (R,t) = - I d(f)cos{<f>R)sm<j) 



tanh(±/?A(/i(i))) 



1 r , w , ,,tanh(i/3A(/i(t))) 

7T 7 A(ft(t)) 

We can now calculate the equilibrium LN E e ^(t) of 
the two-site density matrix p^{t) of the equilibrium state 
p eq {t) via the prescription in Eq. ©. For the same values 
of the parameters as above, so that (3 ps 3.9 via Eq. J2J), 
the equilibrium LN of the two-site equilibrium state has 
an initial value of ps 0.132635. Then E^(t) jumps up 
to w 0.157188 for all t > 0. Note that the equilibrium 
magnetization M z q (t) jumps down in the same situation. 
So with more entanglement, we have more local disorder 
(cf. 0). 

Our interest is to compare this equilibrium two-site 
entanglement with the two-site entanglement of the 
evolved state. The nearest neighbor correlations of the 
evolved state are given by T xy = S(l,t)/i, T xx {t) — 
-G(-M), Tyy[t) = -G(l,t), T zz {t) = A[M z {t)] 2 - 
G(l,t)G{-l,t) + S(l,t)S{-l,t), where G(R,t) and 
S(R,i), for R = ±1, are given by 



G{R,t) 



x I7 2 sin 2 1 



x [{7 2 sin 2 4> 
S(R,t) 



7 r , , tanh(i/3A(a)) 

7T J A(a)A 2 (6) 

(cos (f) — a) (cos (j> — b) 

+(a - fo)(cos </>-&) cos(2A(6)t)] 

1 r , tanh(±/3A(a)) 

7T 7 A(a)A 2 (6) 

(cos 4 1 — a) (cos </> — 6) } (cos <p — b) 

-(a-&)7 2 sin 2 0cos(2A(6)i)], (8) 

7(0-6)1 [* . . ^ sin(2tA(6)) 

#sm(^)sm^ A(fl)A(6) . 



(9) 



The LN Ejy(t) of the two-site density matrix pi2(i) 
of the evolved state p(t) can now be calculated by using 
Eq. J3J. The behavior of LN of the evolved state for 
the same parameters as above, is shown in Fig. ^ With 
time, the entanglement converges to a fixed value. As 7 
approaches 1, the time taken for convergence is longer, 
until for 7=1, the revivals occur for all times. Hence- 
forth, we consider 7 / 1. 

As we see, the two-site entanglement of the state ob- 
tained from the time evolution (from an initial equilib- 
rium state) is not approaching to the two-site entangle- 
ment of the equilibrium state for large times. This sug- 
gests that, just like the (single-site) magnetization, the 
two-site entanglement of an infinite one-dimensional spin 
system in the XY model is nonergodic. 

We will now compare these two nonergodic quantities 
of the XY chain. We will show that their temporal be- 
haviors are in a sense complementary. First note that the 
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FIG. 1: Comparison of temporal behavior of the nearest 
neighbor entanglement to the magnetization, with time, of 
the evolved state of the infinite spin chain in XY model: We 
plot E N (t) and M z (t) for 7 = 0.5, a = 200, a = 0.5, 6 = 0. 



magnetization M z (t) of the evolved state is a damped 
oscillatory function. For our choice of parameters, the 
damping decreases the amplitude of the oscillation, but 
the mean value of the oscillation is more or less fixed 
(see Fig. Therefore, long time average of M z (t) can 
essentially be considered to be a constant at « 0.02, af- 
ter starting with an initial higher value of ps 0.148328. 
Similar feature was seen before for M z q (t). 

The opposite is true when we compare the two-site 
entanglements of the equilibrium state and the evolved 
states. The two-site LN E N (t) of the evolved state con- 
verges to ps 0.18, after starting from an initial lower value 
of pa 0.132635 (see Fig. 0. Again a similar behavior was 
seen in the equilibrium state. It is in this sense that we 
say that the two-site entanglement and (single-site) mag- 
netization has a complementary behavior in the infinite 
spin chain in the XY model. 

The magnetization and nearest neighbor entanglement 
in the evolved state are plotted in Fig. ^ We note here 
that the next-nearest neighbor correlation functions are 
also known [l5T |. Using them, we have calculated the 
LN for the next-nearest neighbor of the evolved state. 
It turns out that it is smaller than the nearest neighbor 
LN. For example, for the parameters as above, the next- 
nearest neighbor LN is vanishing already at t = 0.8 (and 
never going over ps 0.0071), while the nearest neighbor 
LN is about 0.18, on average, for large times. 

Let us now consider the the behavior of the magneti- 
zation and nearest neighbor entanglement of the evolved 
state, for different values of the temperature, but for a 
fixed time t = 1. (See Fig. (2J) As expected, for high tem- 
peratures, the (nearest neighbor) entanglement is vanish- 
ing. Note, that both magnetization and nearest neighbor 
entanglement ultimately saturate, with decreasing T. 

Summarizing, we have presented evidence that in the 
thermodynamic limit of an infinite one-dimensional chain 
of spin- 1/2 particles described by the XY model Hamil- 
tonian, the two-site entanglement is nonergodic. The 
two-site entanglement of the evolved state does not ap- 
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FIG. 2: The behavior in the evolved state, with temperature 
T, of nearest neighbor entanglement and magnetization of the 
infinite spin chain in XY model are compared for a fixed time 
t — 1 and for 7 = 0.5, a = 10, 6 = 0. We actually plot the 
curves with respect to /3, where j3 = 1/fcT. 



proach its equilibrium value. That is, entanglement in 
such systems does not, by itself (i.e. without contact with 
external reservoirs), relax to its equilibrium value, after 
a change of the external magnetic field. This indicates 
that entanglement in such systems cannot be described 
by equilibrium statistical mechanics. We also show that 
the entanglement has a complementary temporal behav- 
ior with respect to magnetization. We believe that such 
studies of the dynamics of entanglement in spin systems 
will help us to implement quantum information process- 
ing tasks in such systems (cf. 8]). 
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